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ABSTRACT 


Work  is  continuing  in  both  gas  phase  and  liquid  phase  turbulent 
mixing  and  chemical  reactions.  In  the  gas  phase  shear  layer  work,  our 
investigations  are  concentrating  on  shear  layer  free  stream  density 
ratio  effects,  finite  kinetic  rate  (Damkohler  number)  effects,  and  heat 
release  effects  the  latter  up  to  a  temperature  rise  in  the  combustion 
zone  of  the  order  of  1  ,000  K.  In  jet  flows,  the  development  of  laser 
Rayleigh  scattering  techniques  is  progressing  towards  conserved  scalar 
measurements  down  to  diffusion  space  and  time  scales.  In  the  liquid 
phase  work,  laser  induced  fluorescence  measurements  in  both  shear  layer 
and  jet  flows  have  yielded  considerable  new  information  on  the  mixing 
process  and  statistics,  permitting  direct  estimates  as  the  probability 
density  function  in  both  cases.  Theoretical  work  in  progress  is 
addressing  the  finite  chemical  rate  problem  as  well  as  the 
diffusion-limited  shear  layer  mixing  problem. 
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1 


1.0  INTRODUCTION 


Work  is  continuing  in  both  gas  phase  and  liquid  phase  turbulent 
mixing,  chemical  reactions  and  combustion.  Both  experimental  and 
analytical/modeling  investigations  are  in  progress,  as  well  as  an 
accompanying  effort  to  develop  and  augment  diagnostic  techniques,  on  an 
as  needed  basis,  in  support  of  the  experimental  program. 


2.0  HF  COMBUSTING  SHEAR  LAYER 

2.1  Density  Ratio  Effects 

Thee  free  stream  density  ratio  across  the  2-D  mixing  layer  has  one 
of  the  largest  effects  on  the  entrainment  and  growth  of  the  turbulent 
region  (Dimotakis  1 98-4 ) .  Using  readily  available  gases,  changes  in  the 
growth  of  the  layer  by  factors  of  three  have  been  previously  reported 
(Konrad  1976).  These  results  have  been  confirmed  in  a  few  pilot  runs  at 
low  heat  release  (Mungal  &  Dimotakis  1984),  but  as  yet  no  systematic 
investigation  of  this  phenomenon  has  been  undertaken. 

In  particular,  we  are  using  the  HF  shear  layer  combustion  facility 
to  determine  how  density  differences  effect  chemical  reactions,  and  in 
the  case  of  high  heat  release,  what  effect  the  existence  of  three 
characteristic  densities  (high  speed  stream,  low  speed  stream,  mixed 
fluid  in  the  layer)  has  on  both  product  formation  and  layer  growth.  The 
facility  was  recently  recalibrated  for  the  changes  in  diluent  gases. 


2.2  Finite  Kinetic  Effects 

An  analysis  of  the  effects  of  finite  chemical  rates  has  been 
carried  out  for  a  number  of  simplified  flow  geometries.  Through  this 
analysis  and  a  model  decomposition  of  the  reacting  zones  within 
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turbulent  shear  layer,  such  as  the  Broadwell-Breidenthal  model,  we  hope 
to  understand  the  data  already  reported,  and  perhaps  predict  these 
effects  for  flows  not  yet  investigated. 

The  simplest  flow  model  analysed  was  the  constant  entrainment, 
perfectly  mixed  region.  For  this  initial  value  problem,  a  simple 
one-step,  second  order  reaction  was  assumed  and  the  temperature 
dependance  of  the  rate  coefficient  was  neglected.  From  this,  both  an 
exact  solution  and  a  number  of  special  case  and  asymptotic  results  have 
been  obtained. 

A  numerical  study  of  the  effects  of  finite  kinetics  on  a  laminar 
strained  diffusion  flame  has  also  been  performed.  In  this  work,  the 
effects  of  equivalence  ratio,  reaction  rate,  strain  rate  and  activation 
energy  on  product  distributions  have  been  investigated.  In  the 
investigation  of  high  activation  energy  effects,  the  hysteresis  curve 
associated  with  the  blow-out  and  light-up  phenomenon  has  been  explored, 
and  the  stability  of  the  solutions  on  this  curve  has  been  calculated. 

This  work  is  carried  out  in  conjunction  with  the  research  described 
in  section  2.5  below. 


2.3  High  Resolution  Temperature  Sensing  Probes 

As  previously  reported,  plans  to  use  the  0.6  ym  temperature  sensing 
probes  have  been  made  and  prototypes  have  been  fabricated  and  tested. 
The  probe  mount,  support  structure  and  purge/ connect  ion  block  have  all 
been  constructed  and  installed.  These  probes,  however,  have  a  reduced 
ability  to  withstand  higher  dynamic  pressures  and  higher  temperatures 
than  our  standard  2.5  ym  probes.  We  have  decided  to  defer  measurements 
using  these  probes,  to  accommodate  the  higher  temperatures  we  would  like 
to  reach  in  the  near  future.  See  section  2.4  below. 
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2.4  Heat  Release  Effects 

High  heat  release  runs  with  adiabatic  flame  temperature  rise  of  up 
to  940  K  are  in  progress.  The  previously  reported  trend  of  layer 
thinning  is  observed  to  persist.  A  substantial  reduction  in  overall 
entrainment  of  free  stream  fluid  into  the  layer  is  caused  by  heat 
release.  It  appears  that  the  reduction  in  entrainment  is  more  than 
sufficient  to  offset  the  displacement  owing  to  density  decrease  in  the 
layer  by  heat  release. 

An  effort  is  being  made  to  model  the  effects  of  heat  release  on  the 
growth  and  entrainment  of  the  layer  using  similarity  methods,  and  it 
seems  that  the  reported  trends  of  layer  thinning  and  entrainment 
reduction  can  be  at  least  partially  accounted  for  by  a  reduction  in  the 
turbulent  shear  stress  by  heat  release. 

Large  scale  structures  are  seen  to  dominate  the  flow  at  all  levels 
of  heat  release.  High  speed  motion  picture  data  were  taken  over  a  range 
of  temperatures  and  are  being  used  to  provide  statistical  information  on 
vortex  spacing.  Preliminary  results  suggest  a  slight  reduction  in  mean 
vortex  spacing  with  heat  release. 

Further  runs  with  a  favorable  pressure  gradient  were  made  at  an 
adiabatic  flame  temperature  rise  of  up  to  553  K  .  For  these  runs  of 
higher  heat  release  and  pressure  gradient,  the  pressure  gradient  still 
does  not  appear  to  have  any  significant  effect  on  the  amount  of  mixing 
and  chemical  product  production  in  the  layer. 

A  more  extensive  discussion  of  these  and  other  results  is  presented 
in  the  appended  AIAA  Paper  85-0142  by  J.  C.  Hermanson,  M.  G.  Mungal 
and  P.  E.  Dimotakis,  recently  presented  at  the  AIAA  23rd  Aerospace 
Sciences  Meeting,  at  Reno,  Nevada.  This  part  of  the  work  is  nearing 
completion.  A  full  documentation  should  be  available  at  the  end  of  the 
next  reporting  period  (ending  April  1986)  in  the  form  of  a  Caltech  Ph.D. 
thesis  by  J.  Hermanson. 
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2.5  Damkohler  number  2  finite  rate  kinetic  effects 

The  chemical  kinetics  of  the  hydrogen/fluorine/nitric  oxide  system 
has  been  studied  using  the  SANDIA  (Livermore)  CHEMKIN  chemical  kinetics 
program.  The  reactions  considered  were 
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Our  nominal  concentrations  investigated  were 

h2  -  4X, 

^2  “0.5/1 

and 

NO  -  0.015$ 

carried  in  a  nitrogen  diluent.  The  overall  chemical  reaction  time  for 
this  system  is  found  to  be  0.5msec. 

At  our  measuring  station,  the  large  scale  mixing  time  is  5msec, 
based  on  visual  width  of  the  shear  layer  and  the  velocity  difference. 
Hence  the  Damkohler  number  (ratio  of  mixing  time  to  chemical  time)  is  10 
when  the  reaction  is  mixing  limited.  This  value  compares  well  with  that 
reported  earlier  by  Wallace  (1981).  The  results  also  show  that  the 
overall  chemical  time  is  controlled  by  the  NO  concentration. 

These  results  are  currently  being  used  to  interpret  the  results  of 
our  attempts  to  simulate  the  effects  of  slow  chemistry.  A  full 


AFOSR-83-0213  Interim  Report 


documentation  (Broadwell  &  Mungal)  on  this  part  of  the  work  should  be 
available  by  the  end  of  the  next  reporting  period.  The  effects  of  slow 
chemistry  are  also  important  in  the  context  of  our  long  term  goals  in 
this  program  as  we  would  hope  to  extend  the  work  on  chemically  reacting 
flows  to  supersonic  flow  speeds,  where  finite  kinetics  play  an  (if  not 
the  most)  important  role. 


3.0  LIQUID  PHASE  SHEAR  LAYER 

Our  laser  induced  fluorescence  (LIF)  measurements  in  the  water  2-D 
mixing  layer  have  shown  that  the  initial  roll-up  of  the  vortical 
structures  is  not  symmetric  in  the  transverse  dimension  of  the  layer. 
Results,  reported  in  our  earlier  progress  reports,  Koochesf ahani  P.h.  D. 
thesis  1984  and  Koochesf ahani  &  Dimotakis,  AIAA  paper  no.  84-0198, 
indicate  that  in  the  initial  roll-up  process  a  large  amount  of  high 
speed  fluid  is  engulfed  into  the  layer,  which  then  leads  into  a  vortical 
structure  whose  core  is  essentially  full  of  fluid  from  the  high  speed 
stream.  The  significance  of  this  event  is  that  the  composition  of  the 
mixed  fluid  during  the  mixing  transition  is  not  determined  solely  by  the 
entrainment  ratio  of  fluids  from  the  two  free  streams  into  the  layer, 
but  also  by  this  initial  composition  asymmetry  in  the  structure  cores  as 
a  result  of  the  first  roll-up. 

There  is  an  effort  currently  under  way  to  try  to  understand  the 
initial  roll-up  process  via  linear  stability  theory  of  the  spatially 
growing  disturbances.  Our  goal  is  to  calculate  the  geometry  of  the 
streakline  that  peels  off  the  splitter  plate  which  subsequently  rolls  up 
into  a  vortex.  The  initial  set  of  calculations,  using  a  hyperbolic 
tangent  mean  velocity  profile  at  the  splitter  plate  as  a  test  case, 
agrees  with  previous  results  of  Michalke  (1965).  See  figure  1  .  We  are 
now  in  the  process  of  calculating  the  streakline  for  the  actual  case 
realized  in  the  laboratory  by  incorporating  the  wake  of  the  splitter 
plate  into  the  mean  velocity  profile. 
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4 . 0  SHEAR  LAYER  MIXING  &  CHEMICAL  REACTIONS  MODEL 

An  analytical  model  of  mixing  and  chemical  reactions,  at  low  heat 
release,  for  fully  developed,  turbulent  shear  layers  has  been  completed. 
The  model  calculates  the  chemical  product  thickness  ratio  6p/g  (volume 
fraction  of  the  total  width  of  the  shear  layer  occupied  by  chemical 
product),  where  6p  is  the  product  thickness  and  6  is  the  visual 
thickness  of  the  shear  layer.  The  resulting  estimates  are  a  function  of 
the  local  shear  layer  Reynolds  number  Re  =  AU  6/v  ,  Schmidt  number 
Sc  >=  v/D,  and  reactant  stoichiometric  mixture  ratio  <{>.  The  overall 
entrainment  ratio  to  the  shear  layer  is  an  input  to  the  calculation,  and 
is  expressed  as  a  function  of  the  velocity  and  density  ratio  of  the 
shear  layer  free  streams,  using  the  results  of  earlier  calculations 
(Dimotakis  1984).  The  model  is  essentially  without  adjustable 
parameters  and  assumes  that  the  spectrum  of  scales  of  the  turbulence 
contained  within  the  confines  of  the  shear  layer  edges  is  described  by  a 
Kolmogorov  turbulent  cascade. 

The  resulting  theoretical  estimates  for  Ap/g,  based  on  this  model, 
are  in  good  agreement  with  the  low  heat  release  data  in  our  HF 
combustion  facility  (Mungal  &  Dimotakis  1984),  as  well  as  the  liquid 
phase  data  of  Koochesf ahani  (1984)  at  a  comparable  Reynolds  number, 
successfully  accounting  for  the  Schmidt  number  effect  of  a  factor  of 
approximately  2  between  the  gas  phase  and  liquid  phase  data.  The  model 
also  correctly  predicted  the  Reynolds  number  dependence  of  the  chemical 
product  in  the  shear  layer,  which  was  recently  measured  by  Mungal, 
Hermanson  &  Dimotakis  (1984). 
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5.0  TURBULENT  JET  MIXING  &  CHEMICAL  REACTIONS 

5. 1  Gas  Phase  Jet  Mixing  ^  Laser  Rayleigh  Scattering 

During  the  past  reporting  period,  the  test  cell  experiments  were 
completed  and  detailed  work  was  begun  on  a  full  scale  experiment 
utilizing  Rayleigh  scattering  techniques.  The  test  cell,  which  has  also 
been  described  in  our  earlier  reports,  has  been  used  to  confirm  our 
performance  calculations  and  provide  insight  into  some  of  the  dynamic 
characteristics  of  our  Rayleigh  scattering  system.  A  schematic  of  the 
test  cell  has  been  included  in  figure  2. 

The  results  of  the  test  cell  experiments  are  good.  The  output 
voltage  from  the  electronics  is  linear  in  the  number  density  of 
molecules  in  the  focal  volume  as  determined  by  the  test  cell  pressure. 
Typical  data  are  given  in  figure  3-  The  small  scatter  in  the  data  is  a 
consequence  of  the  uncertainty  in  the  pressure  reading  and  the  mild 
drift  in  laser  power.  Because  the  gain  of  the  photomultiplier  tube  is 
only  specified  within  a  factor  of  15  by  the  manufacturer  the  slope  of 
the  line  in  figure  3  could  not  be  accurately  predicted.  The 
signal-to-noise  ratio,  however,  can  be  predicted  since  it  is  a  ratio  of 
amplified  quantities.  Since,  in  a  shot  noise  limited  measurement,  the 
noise  voltage  increases  as  the  square  root  of  the  signal  voltage,  the 
noise  voltage  divided  by  the  square  root  of  the  signal  voltage  should  be 
a  constant.  This  "system"  constant,  which  is  related  to  the 
signal-to-noise  ratio,  can  be  predicted  and  measured.  Since  low  values 
of  this  constant  are  optimum,  figure  4  makes  it  clear  that  the 
predictions  fall  about  twelve  percent  below  the  measurements.  Which  is 
reasonable  agreement  since  at  least  3  or  4  multiplicative  factors  in  the 
performance  calculation  can  be  only  be  determined  within  ten  percent. 

After  developing  some  simple  data  acquisition  software  we  made  some 
dynamic  concentration  measurements  in  the  test  cell.  A  transient  jet  of 
high  cross  section  gas  was  introduced  through  a  central  orifice  above 
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the  focal  volume  in  the  test  cell.  The  concentration  history  was 
recorded  using  a  portable  computer.  The  results  showed  that  dust 
particles,  whose  Mie  scattering  cross  sections  are  orders  of  magnitude 
larger  than  the  molecular  Rayleigh  cross  sections  which  comprise  the 
signal,  can  be  readily  recognized  and  intelligent  signal  processing 
should  minimize  their  effect  on  the  concentration  measurements. 

The  design  of  a  full  scale  experiment  has  begun  based  on  the 
promising  results  of  the  test  cell  work.  Laboratory  space  which  has  the 
proper  ventilation,  gas  exhausting  capability,  and  electrical  power 
system  is  being  sought.  Some  machining  has  begun  of  special  optical  and 
mechanical  components. 

This  part  of  the  work  is  co-sponsored  by  the  Gas  Research 
Institute. 

5.2  Liquid  Phase  Jet  Mixing  &  Chemical  Reactions 

Our  work  on  entrainment,  mixing  and  chemical  reaction  in  high 
Schmidt  number  (liquid  phase)  turbulent  jets  has  continued. 

We  have  completed  direct  measurements  of  concentration  using  LIF 
digital  imaging  simultaneously  at  1024  points  along  the  jet  axis.  These 
measurements  have  yielded  axial  PDFs  of  jet  fluid  concentration,  and 
have  revealed  a  similarity  concentration  in  terms  of  which  the  PDFs  in 
the  jet  display  self-similarity  along  any  ray  (  r/x  =  constant  ) .  This 
self-similarity  permits  the  cross-over  from  neat — field  to  far-field  jet 
mixing  to  be  defined  quantitatively,  and  allows  the  self-similar  form  of 
the  PDF  throughout  the  jet  far-field  to  be  determined  without  the 
resolution  limitations  inherent  in  any  direct  measurement. 

We  have  also  completed  similar  imaging  measurements  of 
instantaneous  radial  concentration  profiles  across  the  jet  in  the  far 
field  (see  Dahm  &  Dimotakis  1985,  appended).  These  results  demonstrate 
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quantitatively  that  instantaneous  profiles  of  concentration  do  not 
resemble  the  mean  concentration  profile.  Specifically,  the  mixed  fluid 
within  large  regions  in  the  jet  usually  lies  within  a  narrow  range  of 
values,  and  unmixed  ambient  fluid  can  be  found  throughout  the  jet.  The 
mean  profile  is,  therefore,  a  poor  representation  of  the  mixed  fluid 
states  within  the  jet.  This  contrasts  with  the  classical  picture  of  jet 
mixing,  and  may  have  particularly  important  implications  for  the 
modeling  of  jet  mixing  and  the  prediction  of  chemical  reactions  between 
the  mixed  and  ambient  entrained  fluid. 

A  full  documentation  of  this  effort  should  be  available  in  the  next 
reporting  period  (ending  15~Apr-86)  in  the  form  of  a  Caltech  Ph.D. 
thesis  by  W.  Dahm. 

This  part  of  the  work  is  co-sponsored  by  the  Gas  Research 
Institute. 

5.3  Laser  Soot  Scattering  in  Jets  and  Buoyant  Flames 

A  novel  combustion  diagnostic  technique  is  being  explored  which 
makes  use  of  laser  light  scattering  from  soot  particles  in  a  flame. 
While  many  laser-based  diagnostics  are  limited  by  the  presence  of  soot, 
the  soot  particles  themselves  can  in  fact  be  utilized  as  an  important 
marker  within  the  flame.  The  submicron  particles  are  strong  scatterers 
of  visible  radiation  (relative  to  molecular  scattering  cross-sections) 
and  by  directing  the  intense,  monochromatic  beam  of  a  laser  through  the 
soot  distribution  (flame),  the  location  of  soot  along  the  line  of  the 
beam  can  be  determined.  The  essence  of  this  new  technique  is  to  measure 
the  intensity  of  the  soot  scattering  along  a  specified  line  segment  in 
the  flame  as  a  function  of  time. 

Photodiode  array  and  computer  data  acquisition  technology  for 
measuring  and  recording  such  an  intensity  pattern  has  been  developed  and 
extensively  used  in  our  laboratories.  Its  application  in  a  soot 
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scattering  experiment  is  straightforward.  Some  preliminary  feasibility 
studies  using  a  3  watt  argon-ion  laser  and  a  buoyant,  sooting  flame 
indicated  that,  by  visual  inspection,  measurable  scattered  light  appears 
to  be  present.  Although  flame  luminosity  could  be  expected  to  interfere 
with  such  a  measurement,  the  very  narrow  line  width  of  the  laser 

combined  with  a  narrow  band  optical  filter  greatly  diminishes  this 
problem.  This  technique  should  contribute  to  a  variety  of  combustion 
studies,  potentially  addressing  such  questions  as  flame  structure, 

combustion  zone  location  and  character,  and  the  process  of  sooting. 

In  order  to  make  more  quantitative  measurements  and  to  get  some 
idea  of  the  type  of  information  that  can  be  gained  with  such  a 
technique,  an  optical  system  was  designed.  The  purpose  of  the  optical 
system  is  to  pass  the  laser  beam  through  the  flame  in  an  easily 

controlled  manner  and  to  collect  only  the  scattered  laser  light, 
discriminating  against  the  unwanted  flame  luminosity.  With  this 
apparatus,  several  sets  of  photographs  of  the  buoyant  methane  flame 
above  a  50  cm  diameter  glass-bead  burner  were  taken.  These  photographs 
verify  our  initial  estimates  of  the  technique's  efficacy  and  provide  a 
few  insights  into  the  soot  distribution  patterns  within  the  flame. 

These  photographs  indicate  that  scattered  radiation  is  always 

closely  associated  with  a  specific  intersection  of  the  luminous  soot 
distribution  with  the  laser  beam.  We  saw  no  evidence  of  appreciable 
laser  scattering  from  any  region  where  there  was  no  luminosity.  In 
other  words,  we  saw  no  appreciable  cold  soot  far  from  a  flame  sheet.  In 
this  sense,  the  soot  can  be  considered  as  approximately  marking  parts  of 
the  flame  sheet  surface. 

Perhaps  even  more  striking  was  the  observation  that  the  soot  fills 
a  very  small  fraction  of  the  total  volume  of  the  flame.  When  the  soot 
distribution  intersects  the  laser  beam  segment,  it  showed  up  as  a  few 
small  points  along  the  segment,  occasionally  appearing  as  a  short  length 
much  shorter  that  the  width  of  the  overall  flame.  These  observations 
are  consistent  with  a  very  convoluted  flame  sheet  that  has  been  strained 
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out  (at  least  as  far  as  soot  production  is  concerned)  at  some  regions  in 
its  topology  and  which  intersects  the  segment  under  study  at  only  a  few 
points  or  occasionally  as  a  short  segment  when  the  angle  of  intersection 
is  quite  oblique. 

We  are  presently  working  on  implementing  the  photodiode  array  to 
image  the  segment,  using  the  same  beam  steering  and  collection  optics. 
By  interfacing  the  array  to  a  computer  data  acquisition  system  with  high 
speed  A/D  conversion,  we  will  record  the  soot  distribution  along  the 
imaged  segment  as  a  function  of  time.  This  will  give  us  a  more  complete 
picture  of  the  structure  of  the  soot  distribution  and  the  flame  sheet. 

This  part  of  the  work  is  co<-sponsored  by  the  Gas  Research 


Institute. 
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ABSTRACT 

The  effects  of  heat  release  are  studied  In  a 
planar,  gaseous  reacting  mixing  layer  formed 
between  free  streams  containing  hydrogen  and 
fluorine  in  inert  diluents.  Sufficiently  high 
concentrations  of  reactants  are  employed  to 
produce  adiabatic  flame  temperature  rises  of  up  to 
990  K  (1290  K  absolute).  Though  the  displacement 
thickness  of  the  layer,  for  zero  streamwlse 
pressure  gradient,  increases  with  increasing  heat 
release,  the  actual  thickness  of  the  mixing  layer 
is  not  observed  to  Increase,  and  is  in  fact 
characterized  by  a  slight  thinning.  The  overall 
entrainment  into  the  layer  is  seen  to  be 
substantially  reduced  by  heat  release.  The  large 
scale  vortical  nature  of  the  flow  appears  to 
persist  at  high  temperatures.  Imposition  of  a 
favorable  pressure  gradient,  though  resulting  in 
additional  thinning  of  the  layer,  is  observed  to 
have  no  resolvable  effect  on  the  mixing  and  amount 
of  chemical  product  formation. 


INTRODUCTION 

This  investigation  is  concerned  with  heat 
release  effects  in  a  gas  phase  turbulent,  plane, 
reacting  shear  layer  at  high  Reynolds  number.  The 
work  is  an  extension  of  earlier  work  in  the  same 
f  acili  ty 1  ^ -20 ,21  .  The  now  consists  of  a 
two-dimensional  mixing  layer  with  gas  phase  free 
streams,  one  stream  carrying  a  given  concentration 
of  hydrogen  in  an  inert  diluent;  the  other, 
fluorine.  The  reaction  H2  *  f2  *  2  HF  is  highly 
exothermic  30  that  reactant  concentrations  of  11 
h2  and  11  F2  in  an  N2  diluent  respectively  produce 
an  adiabatic  flame  temperature  rise  of  93  K  above 
ambient.  Results  will  be  presented  here 
corresponding  to  fluorine  concentrations  of  up  to 
61  and  hydrogen  concentrat tons  of  up  to  2M1  ,  with 
a  maximum  adiabatic  flame  temperature  rise  of 
990  K  (corresponding  to  1290  K  absolute). 

In  earlier  work  by  Mungal  et  al'^'21^2',  with 
temperature  rises  of  up  to  165  K,  no  coupling  of 
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heat  release  with  the  fluid  mechanics  could  be 
observed,  as  manifested  by  the  growth  rate, 
entrainment  and  discernible  large  3cale  structure 
dynamics.  In  that  work,  the  chemical  reaction 
could  be  considered  a3  a  diagnostic  to  infer  the 
amount  of  molecular  mixing  without  disturbing  the 
overall  properties  of  the  layer.  In  the  work 
reported  here,  the  heat  release  is  much  larger  and 
the  effects  of  the  heat  release  itself  on  the 
properties  of  the  shear  layer  are  investigated. 
The  highest  heat  release  cases  reported  here  are 
sufficiently  hot  to  result  in  a  mean  density 
reduction  in  the  center  of  the  layer  of  about 
65*  . 


EXPERIMENTAL  FACILITY  AND  INSTRUMENTATION 

The  experimental  apparatus  is  described  in 
detail  in  Mungal  et  al20.  It  i3  a  blowdown 
facility  in  which  prenuxed  volumes  of  fluorine  in 
an  inert  diluent  and  hydrogen  in  an  inert  diluent 
are  discharged  through  sonic  orifices  to  maintain 
a  constant  ma33  flux.  Each  stream  enters  a 
settling  and  contraction  section  for  turbulence 
suppression  with  the  high  speed  stream  emerging 
from  a  6:1  contraction  with  an  exit  area  of 
5  «  20  era,  and  the  low  speed  emerging  from  a  U;1 
contraction  in  a  7.5  « 20  cm  exit  area.  The  two 
streams  meet  at  the  tip  of  a  splitter  plate,  with 
a  trailing  edge  Included  angle  of  3.78*.  The  nigh 
speed  free  stream  turbulence  level  was  measured  to 
be  about  2/3*  . 

To  offset  the  free  stream  density  difference 
that  results  from  large  amounts  of  hydrogen  in  one 
stream,  the  densities  of  the  free  streams  were 
matched  for  most  cases  by  using  as  diluent  a 
mixture  of  nitrogen  and  a  small  amount  of  helium, 
on  the  fluorine  side,  and  a  mixture  of  nitrogen 
with  a  small  amount  or  argon  on  the  hydrogen  side. 
The  Intent  was  also  to  match  the  heat  capacities 
of  the  two  free  streams. 

Runs  were  performed  with  a  nominal  high  speed 
flow  velocity  of  22  m/s  ,  and  a  free  stream  speed 
ratio  of  l^/U-]  *0.9  .  In  practice,  these  values 
varied  slightly  —  a  result  of  differences  in  gas 
constants  of  the  various  mixtures  --  although  the 
sonic  metering  orifices  were  adjusted  for  each  run 
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to  minimize  these  variations.  The  measuring 
station  was  positioned  45.7  era  downstream  of  the 
splitter  plate  trailing  edge.  The  Reynolds  number 
at  the  measuring  station  was  typically 
Re5  -6  *10y»  based  on  the  free  stream  velocity 
difference,  the  !X  thickness  of  the  mean 
temperature  profile,  and  the  cold  free  stream 
kinematic  viscosity.  The  IX  thickness,  61  ,  of 
the  temperature  or  concentration  field  1®»20,  13 
defined  here  as  the  transverse  width  of  the  layer 
at  which  the  mean  temperature  rise  is  IJ  of  the 
maximum  mean  temperature  rise  and  has  been  shown2*-* 
to  correlate  well  with  the  visual  thickness  <5Vis 
(Brown  &  Roshko^)  of  the  layer.  Thi3  value  of  the 
Reynolds  number  is  well  above  that  for  the  mixing 
transition  as  reported  by  Bernal  et  al^, 
Breidenthal^  and  Konrad1*^.  The  corresponding 
Reynolds  number  based  on  the  high  speed  free 
stream  velocity  and  on  the  downstream  distance  was 
Rex«6«10^.  A  diagram  of  the  shear  layer 
geometry  is  shown  in  Figure  1 . 


Fig.  I  Turbulent  Shear  Layer  Geometry 


Temperature  data  were  recorded  with  a  rake  of 
eight  2.5  urn  diameter  plat inum  -  1  OX  rhodium  cold 
wires,  with  a  typical  wire  span  of  1.5mm,  welded 
to  Inconel  prongs.  For  some  runs  a  rake  of  25  um 
dlimeter  Chromel-Alumel  thermocouples  was 
employed.  It  was  found  that  the  2.5  um  resistance 
wires  in  the  hottest  regions  did  not  survive  in 
runs  in  which  the  adiabatic  flame  temperature  rise 
exceeded  approximately  600  K.  Both  the  cold  wire 
and  thermocouple  rakes  were  positioned  across  the 
transverse  extent  of  the  layer.  The  total  data 
rate  for  the  resistance  wires  was  30  kHz, 
correspond  1 ng  to  10  kHz  per  probe.  The 
thermocouples  were  sampled  at  500  Hz  each,  for  a 


total  data  rate  of  4  kHz;  their  considerably 
lower  frequency  response  not  warranting  a  higher 
rate . 

Thermocouples  produce  a  voltage  proportional 
to  the  Junction  temperature  and  normally  do  not 
require  calibration.  The  resistance  wires  were 
calibrated  a3  described  previously1^  using  a  hot 
and  cold  Jet  of  known  temperature.  The  two 
measurements  provided  calibration  constants  to 
convert  voltage  to  temperature  rise.  An 
additional  correction  was  applied  to  the  output 
signal  voltage  in  the  present  experiments  to 
account  for  the  non-linearity  in  the  resistivity 
of  the  platinum  -  1 0X  rhodium  wire  element  at 
elevated  temperatures1 ^ .  It  was  determined  that 
for  neither  the  thermocouples  nor  the  resistance 
wires  was  there  significant  radiation  error  for 
the  temperatures  in  this  Investigation214.  Both 
probes,  however,  are  influenced  by  heat  conduction 
to  the  support  prongs,  which  can  result  In 
excursions  from  the  mean  temperature  being  in 
error  by  as  much  as  10X  to  20X  for  the  cold  wires 
and  up  to  40X  for  the  thermocouples.  Both 
diagnostics,  however,  produce  accurate  mean 
temperatures,  as  during  a  small  fraction  of  the 
course  of  the  run  (before  data  acquisition  begins) 
the  tips  of  the  support  prongs  assume  the  local 
mean  value.  Good  agreement  (typically  within  5X) 
was  obtained  during  runs  In  which  both  sets  of 
probes  were  employed.  Errors  resulting  from 
differences  in  the  thermal  conductivity  of  the 
fluid  from  each  of  the  free  streams  were 
established  to  be  small. 

In  addition  to  the  temperature  data,  a 
schlieren  system  was  utilized  for  concurrent  flow 
visualization.  The  beam  width  utilized  was 
sufficient  to  illuminate  25  cm  of  the  shear  layer. 
A  circular  source  mask  and  a  circular  hole  spatial 
filter  were  used  in  place  of  the  conventional 
source  3llt  and  knife  edge  In  an  effort  to  give 
equal  weights  to  gradients  in  index  of  refraction 
in  all  directions  and  better  resolve  the  large 
scale  structure  of  the  flow.  The  hole  sizes  were 
increased  with  increasing  flow  temperature  to 
optimize  (reduce)  sensitivity  as  needed.  High 
time  resolution  spark  schlieren  photographs  were 
taken  with  a  spark  source  (-  3  usee  duration), 
synchronized  with  a  motor-driven  35  mm  camera,  at 
a  rate  of  approximately  three  frames  per  second. 

The  mean  velocity  profile  was  measured  for 
most  run3  by  a  pitot  probe  rake  of  fifteen  probes 
connected  to  a  miniature  manometer  bank  filled 
with  fluorine  resistant  oil  (Hooker  Chemical 


Fluorolube  FS-5)  with  an  adequate  time  response  to 
yield  a  reliable  mean  dynamic  pressure  profile 
during  each  run.  The  bank  was  photographed  by  a 
second  motor-driven  35  mm  camera.  The 
photographic  data  were  digitized  and  reduced  to 
mean  velocity  profiles.  This  technique  of 
measuring  the  pitot  pressure  was  estimated  to  be 
accurate  to  5t.  Rebollo2^  estimated  that  the 
accuracy  of  extracting  mean  velocities  from  pitot 
pressures  in  non-constant  density  flows  is  about 
9-5  t  . 

Finally,  the  3treamwise  static  pressure 
gradient  was  monitored  by  measuring  the  pressure 
difference  between  two  downstream  locations  on  the 
low  speed  sidewall  with  a  Datametrlcs  type  573 
fluorine  resistant  Barocel  sensor.  The  high  speed 
side  wall  was  kept  horizontal  for  all  runs  and  the 
low  speed  side  wall  was  adjusted  for  the  desired 
streamwlse  pressure  gradient.  The  wedge-like 
geometry  of  the  planar  shear  layer  displacement 
allows  this  simple  means  of  accommodating  or 
imposing  any  desired  pressure  gradient.  Most  of 
the  runs  in  the  present  investigation  were 
performed  with  the  sidewalls  adjusted  to  the 

requisite  divergence  angle  to  ensure  a  zero 

streamwlse  pressure  gradient.  For  some  runs  at 

high  heat  release,  the  walls  were  left  fixed  at 

the  angle  required  for  zero  pressure  gradient  at 
zero  heat  release,  which  allowed  the  heat  release 
to  induce  a  favorable  streamwlse  pressure  gradient 
(accelerating  flow),  as  a  result  of  the  combustion 
displacement  effects. 

CHEMISTRY 

The  chemical  reaction  utilized  in  the  present 
investigation  is  effectively 

♦  F2  -  2HF  ,  AQ  -  -  130  kcal/mole 

Thi3  yields  a  temperature  rise  of  93 K  for  1%  ?2 
and  11  H,  in  N2  diluent  under  constant  pressure, 
adiabatic  conditions  (this  is  the  so-called 
adiabatic  flame  temperature  rise).  The  chemical 
reaction  actually  consists  of  two  second  order 
chain  reactions: 

Hj  »  F  ♦  HF  *  H,  a Q  *  -  32  kcal/mole 

H  *  *  HF  *  F,  AQ  *  -  98  kcal/mole 


Proper  chain  initiation  requires  3ome  free  F 
atoms,  which  were  generated  in  these  experiments 
by  premixing  a  trace  amount  of  nitric  oxide  into 
the  hydrogen-carrying  stream.  This  allows  the 
reaction 

F2  *  NO  »  NOF  *  F 

which  provides  the  required  small  F  atom 
concentration  in  the  layer  to  sustain  proper 
ignition.  For  all  runs  in  this  investigation,  the 
NO  concentration  was  maintained  at  31  of  the  free 
stream  fluorine  concentration. 

For  all  flows  reported  here,  the  resulting 
chemical  time  scales  are  fast  compared  with  the 
fluid  mechanical  time  scales.  The  chemical  time 
scales  for  the  reaction,  over  the  entire  range  of 
concentrations ,  were  determined  using  the  CHEMKIN 
chemical  kinetics  program1 5.  The  chemical  rate 
data  for  the  reactions  involved  were  taken  from 
Cohen  !>  Bott11  and  Baulch  et  al2.  The  Dahnkohler 
number  (ratio  of  mixing  time  to  chemical  time) 
based  on  the  local,  large  scale  characteristic 
time  (A1  /AU,  where  AU  is  the  free  stream 
velocity  difference)  ranges  from  25  to  125  with 
increasing  reactant  concentrations.  The 
Dahmkohler  number  based  on  the  time  of  flight  from 
the  splitter  tip  to  the  measuring  station  (  x/Ug, 
where  Uc  is  the  large  structure  convection 
velocity)  varies  from  l 60  to  800.  For  comparison, 
the  Dahmkohler  numbers  for  the  low  heat  release 
work  of  Mungal1^’20’21  were  10  and  60  for  the 
large  scales  and  time  of  flight,  respectively.  In 
those  investigations,  the  chemistry  was  shown  to 
be  already  fast.  Chemical  kinetics  are, 
consequently,  not  an  issue  in  the  present 
investigation,  where  the  chemistry  is  much  faster 
as  a  result  of  the  higher  reactant  concentrations 
and  combustion  temperatures. 

The  equivalence  ratio,  *,  is  defined  here  as 
the  volume  of  high  speed  fluid  required  to 
completely  react  with  a  unit  volume  of  low  speed 
fluid.  This  is  the  same  as  the  ratio  of  the  low 
speed  free  stream  molar  concentration,  Cq2,  to  the 
high  speed  free  stream  molar  concentration,  c^, , 
divided  by  the  low  speed  to  high  speed 
stoichiometric  ratio,  l.e. 

;02/c01 


In  thl3  case,  since  the  molar  stoichiometric  ratio 
for  the  hydrogen-f luor ine  reaction  is  unity. 
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RESULTS  AND  DISCUSSION 


I .  Growth  Rate  and  Entrainment 

The  low  speed  sidewall  divergence  required 
for  zero  streamwlse  pressure  gradient  13  a  direct 
measure  of  the  displacement  thickness  of  the 
layer,  6  ,  where  6*/x  indicates  the  tangent  of 
the  angle  by  which  the  low  speed  free  streamline 
Is  shifted  owing  to  the  presence  of  the  3hear 
layer.  Note  that  the  displacement  thickness  Is 
less  than  zero  for  a  layer  with  no  heat  release, 
and  increases  3teadily  with  heat  release,  as  shown 
in  Figure  2  .  The  parameter  (  pQ  -  p  )/pQ 
represents  the  mean  normalized  density  reduction 
in  the  layer  due  to  heat  release,  where  p  Is  the 
mean  density  in  the  layer  and  pQ  Is  the  average 
(cold)  density  of  the  free  streams. 


with  the  identical  3peed  and  density  ratio,  using 
a  formula  derived  in  Dimotakis12  .  A  straight 
line  least  squares  fit  to  the  data  suggests  that 
the  thinning,  for  a  mean  density  reduction  of 
40t  ,  may  be  as  high  as  15'  .  No  dependence  of  tne 
thinning  trend  on  equivalence  ratio  was  observed. 


Fig.  3  Leyoc  Thinning  Mith  Hoot  Aalesaa 

The  slight  reduction  in  layer  thickneso  with 
increasing  heat  release  13  confirmed  by  the  mean 
velocity  data.  Sample  velocity  profiles,  at 
different  heat  release  but  identical  speed  and 
density  ratios,  are  presented  in  Figure  4  .  it 
can  be  seen  that  the  hotter  layer  is  noticeably 
steeper  in  maximum  slope,  in  agreement  with 


Effect  ol  Maat  Releeee  on 
Shaar  Layer  Maan  Valoclty  Profile 


It  may  be  worth  noting  that  the  actual  shear 
layer  thickness  (6^  or  <5via),  In  spite  of  large 
heat  release  and  large  density  changes,  does  not 
increase  and,  in  fact,  shows  a  slight  decrease, 
even  though  the  displacement  thickness  Increases 
with  heat  release.  This  effect  was  noted  by 
Wallace2^  and  is  observed  in  the  present  set  of 
experiments,  in  which  the  maximum  mean  flow 
temperature  increase  is  about  three  times  greater 
than  that  of  Wallace.  The  observed  1  %  temperature 
profile  thickness  at  zero  pressure  gradient  is 
plotted  versus  the  mean  density  in  the  layer  in 
Figure  3-  Since  it  is  difficult  to  hold  the  speed 
ratio  at  exactly  0.40  from  run  to  run  and  also 
because  the  density  ratio  of  tne  free  streams  is 
slightly  different  from  unity  for  some  runs,  each 
data  point  has  been  corrected  by  normalization 
with  the  expected  growth  rate  for  a  cold  layer 


v/ X 


Effect  ol  Hoot  Releeee  on 
Shaar  Layer  Moon  Velocity  Profile 


Wallace^.  Normalization  of  this  maximum  slope  by 
the  free  stream  velocity  difference  gives  the 
vorticity  thickness  of  the  layer. 


A  plot  of  the  vorticity  thickness  variation  with 
heat  release,  again  corrected  for  variations  in 
speed  ratio  and  density  ratio,  is  shown  in  Figure 
5  .  Each  point  is  normalized  by  the  corresponding 
cold  flow  vorticity  thickness.  The  portion  of  the 
present  results  at  moderate  heat  release, 

including  some  of  the  unpublished  data  of  Mungal, 
are  seen  to  be  in  good  agreement  with  Wallace^. 
The  abscissa  parameter  of  Figure  5*  9max 
ATmax/TQ  ,  the  maximum  mean  temperature  rise  over 
the  ambient  temperature,  scales  differently  than 
(p0-p)/pQ.  the  normalized  mean  density 

reduction,  and  renders  the  thinning  effect  at 
lower  temperatures  more  apparent. 


Xq,  since  the  relevant  similarity  downstream 
coordinate  is  in  fact  y  /  (  x  -xQ)  .  The  trends  in 
layer  thinning  reported  here  do  allow  the 
possibility  that  some  of  the  effects  could  be 
accounted  for  by  a  3hift  in  the  virtual  origin 
with  heat  release.  The  changes  in  xQ  with  heat 
release  have  not  been  investigated  as  yet.  That 
initial  conditions  can  have  a  significant  effect 
on  layer  growth  has  been  shown,  for  example,  by 
Browand  &  Latigo^.  See  also  Batt1,  Bradshaw14  and 
the  discussion  in  the  review  paper  by  Ho  & 
Huerre1 44 . 

One  implication  of  the  fact  that  the  layer 
width  does  not  increase  with  increasing 
temperature  is  that,  since  the  density  in  the 
layer  is  substantially  reduced  but  the  layer  does 
not  grow  faster,  the  volumetric  entrainment  of 
free  stream  fluid  into  the  layer  must  also  be 
greatly  reduced  by  heat  release.  The  amount  of 
entrainment  into  the  layer  can  be  calculated  from 
the  mean  velocity  and  density  (i.e.  temperature) 
profiles  as  follows: 


Fig.  3  Effect  of  nalassa  on  vorticity  TMcfcn««a 


An  Investigation,  at  higher  temperatures,  by 
Pitz  i  Daiiy^2  in  a  combusting  nixing  layer  formed 
downstream  of  a  rearward  facing  step  found  that 
the  vorticity  thickness  did  not  appear  to  change 
between  their  cold  runs  and  hign  heat  release 
runs.  Keller  &  Daily13,  however,  report  that  in  a 
reacting  mixing  layer  between  a  cold  premixed 
reactant  stream  and  a  preheated  combustion  product 
3tream,  the  vorticity  thickness  increases 
signif lcantiy  with  increasing  temperature.  The 
reasons  for  the  discrepancy  between  these  results 
and  the  ones  reported  here  are  not  clear  at  this 
writing. 

A  complicating  factor  in  ar  discussion  of 
growth  fate  is  the  location  of  tne  virtual  origin. 


where  V  Is  the  volume  flux  into  the  layer  per  unit 
span,  x  is  the  downstream  coordinate,  and  n  -  y/x 
is  the  shear  layer  similarity  coordinate.  Thi3 
expression  assumes  that  the  layer  is  self-similar 
at  the  station  at  which  tne  integral  is  performed. 
The  quantity  pU  was  computed  as  p  U,  which  was 

used  here  as  an  approximation  for  the 

density-velocity  correlation  pU. 

Results  'Yom  Mungal  et  al^1  suggest  that 
there  is  a  Reynolds  number  dependence  on  product 

formation.  Since  the  growth  rate  does  appear  to 

be  a  function  of  the  product  formation  (i.e.  heat 
release),  strictly  speaking,  the  flow  would  not  be 
expected  to  be  exactly  self  similar.  An  alternate 
method**  *■  which  approximates  the  overall 

entrainment  is  to  usep  the  geometry  of  tne  layer  as 
3hown  In  Figure  1  to  derive: 

V 

77 —  -  ni  “  r  ( n?  *  tan  8)  , 

Uj  x 

wnere  r  -  »  nt,2  are  similarity 

coordinate  edges  of  the  shear  layer,  and  3  is  the 
deflection  angle  of  the  low  speed  aide  wall.  A 
common  difficulty  of  both  methods  13  that  of 
selecting  the  proper  values  for  m  and  ^  .  One 
reasonable  choice  is  the  pai"  of  points 
corresponding  to  the  1%  edges  of  the  temperature 


profile.  Resulting  calculations  for  choices  of  1% 
and  10%  points  in  the  temperature  profiles,  for 
both  the  Integral  and  geometric  methods,  are 
plotted  in  Figure  6  .  It  can  be  seen  that, 
regardless  of  the  choice  of  edge  reference  points, 
the  Inference  is  that  the  entrainment  into  the 
layer  is  strongly  reduced  as  a  function  of  heat 
release,  amounting  to  about  50%  ,  for  a  mean 
density  in  the  layer  of  40%  below  its  nominal  cold 
value.  That  the  entrainment  reduction  is  in 
exccess  of  the  mean  density  reduction  suggests  that 
the  'decrease  in  entrainment  flux  more  than 
compensates  for  the  additional  displacement  due  to 
density  change. 


IPo-flfPo 


Fig.  0  Effect  of  Moot  Boloooo  on  Overall  Layer  Entralmaent 


II .  Large  Scale  Structure  Dynamics 

Figures  7a  and  7b  show  time  traces  of 
temperature  from  the  rake  of  cold  wires.  Both 
runs  are  at  *  •  I.  The  data  in  Figure  7a  are  Trom 
a  run  with  2%  F2  and  2%  H2,  for  a  flame 
temperature  ri3e  of  186  K,  Figure  7b  shows  the 
results  of  a  run  with  6%  F2  and  6%  H2,  with  a 
flame  temperature  rise  of  553  K.  At  elevated 
temperatures ,  the  flow  dynamics  still  appear  to  be 
dominated  by  large  scale  structures  separated  by 
cold  tongues  of  fluid  which  extend  well  into  the 
layer.  These  findings  are  consistent  with  the 
earlier  results  at  low  temperatures 1 9 > 20, 21 . 


Fig.  7  Teepereture  v«.  Time  Trecee 


Spark  schlieren  photographs  are  presented  in 
Figures  8a  and  8b,  representing  conditions  of  both 
low  and  high  heat  release.  Figure  8a  is  a  2%  F2 
and  2%  H2  run  with  an  adiabatic  flame  temperature 
rise  of  186  K.  I  Figure  8b  is  shown  a  6%  F,  and 
24%  H2  run  with  an  adiabatic  flame  temperature 
rise  of  847  K.  The  large  scale  structure  is 
evident  in  both  photographs  and  does  not  appear  to 
be  greatly  altered  by  the  heat  release. 
Furthermore,  that  the  structures  appear  well 
defined  in  the  schlieren  photographs  suggests  that 
they  retain  their  predominately  two-dimensional 
nature.  Ganj  i  &  Sawyer1^  and  Pitz  &  Daily22 
continued  to  observe  the  large  scale  structures 
behind  a  rearward  facing  step  at  higher 
temperatures  still,  corresponding  to  adiabatic 
flame  temperatures  of  up  to  about  1650 K  absolute. 
Keller  &  Dally10  al30  observed  that  well-ordered 
structures  were  present  for  all  values  of  heat 
release  in  their  investigations. 
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Fig.  8  Spark  Schlieren  Photographs 
of  Reacting  Shear  Layers 


III.  Pressure  Gradient  Results 

It  could  be  argued  that.  In  combusting  flows 
with  substantial  density  variations  and 
appreciable  pressure  gradients.  an  additional 
mixing  mechanism  might  be  operative,  resulting 
from  a  possible  relative  acceleration  between 
light  fluid  elements  and  heavy  fluid  elements^ • . 
The  efficacy  of  such  a  mechanism  would  of  course 
depend  on  the  scales  at  which  the  density 
variations  would  be  observed  and  their  relation  to 
the  viscous  small  scales  of  the  flow.  If  the 
hot/cold  fluid  elements  are  very  closely  spaced, 
viscous  effects  might  not  permit  large  relative 
motions  to  be  established  and  little  or  no 
augmentation  of  the  mixing  would  be  observed. 
These  issues  were  investigated,  in  part,  In  the 
experiments  discussed  below. 

Setting  the  sidewalls  for  zero  pressure 
gradient  for  the  cold  flow  yields  a  naturally 
induced  favorable  pressure  gradient  in  the  case  of 
flow  with  heat  release.  In  thi3  investigation, 
such  favorable  pressure  gradient  runs  were  made  at 


reactant  concentrations  up  to  6»  and  6X  Fj, 
corresponding  to  a  flame  temperature  rise  of 
553  K.  This  was  sufficient  to  induce  a  pressure 
increment  sp,  over  the  distance  from  the  splitter 
plate  tip  to  the  measuring  station,  of  about  1/2 
of  lpQ  Uj  ,  the  low  speed  free  stream  dynamic 
head.  This  caused  the  high  speed  velocity  to 
Increase  from  21.2  to  22.3  m/s  and  the  low  3peed 
velocity  from  8.2  to  10.3  m/s  between  the  splitter 
tip  and  the  measuring  station. 

The  resulting  mean  temperature  profiles,  with 
and  without  pressure  gradient,  are  presented  in 
Figure  9.  The  layer  is  thinned  by  the  pressure 
gradient,  a3  would  be  expected  in  accelerating 
flow.  It  can  be  seen,  however,  at  least  for  the 
values  of  heat  release  and  pressure  gradient 
reported  here,  that  there  appears  to  be  no 
significant  change  in  either  the  peak  temperature 
or  the  total  amount  of  product  formation  Induced 
by  the  favorable  pressure  gradient.  Temperature 
profiles  at  lower  temperatures  also  3how  no 
systematic  changes  resulting  from  pressure 
gradient,  at  least  to  within  the  estimated 
reproducibility  and  accuracy  of  the  data  (3"5X). 


BX  F2  GX  H2  ATf  >  553  K 
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Fig.  9  Effect  of  Pressure  Gradient 
on  Mean  Temperature  Profile 
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CONCLUSIONS 

We  have  found  that  the  growth  rate  of  a 
chemically  reacting  shear  layer  with  heat  release 
does  not  Increase  and  in  fact  appears  to  decrease 
slightly.  In  the  presence  of  an  Increase  in  the 
shear  layer  displacement  thickness  as  a  result  of 
heat  release,  one  might  have  expected  a 
commensurate  increase  in  shear  layer  thickness. 
The  Implication  is  that  the  decrease  in  the 
entrainment  flux  must  more  than  compensate  for  the 
displacement  effect. 

The  imposition  of  a  favorable  pressure 
gradient  13  not  found  to  have  any  noticeable 
effect  on  the  amount  of  mixing  and  chemical 
production  in  the  layer. 
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Abstract 

An  experimental  Investigation  of  entrainment 
and  mixing  In  the  self-similar  far  field  of  an 
axlsymmetrlc  free  turbulent  Jet  In  water  is 
presented.  Length  and  time  scales  for  the  flame 
length  fluctuations  of  reacting  jets  are  shown  to  be 
approximately  equal  to  the  local  characteristic 
large  scale  length  and  time  of  the  flow.  It  Is  also 
shown  that  Instantaneous  radial  profiles  of 
concentration  across  the  Jet  do  not  resemble  the 
mean  concentration  profile,  Indicating  that  the  mean 
profile  Is  a  poor  representation  of  the  mixed  fluid 
states  within  the  Jet.  These  instantaneous  profiles 
also  show  that  unmlxed  ambient  fluid  Is  transported 
throughout  the  entire  extent  of  the  Jet,  and  that 
the  mixed  fluid  composition  within  the  jet  can  be 
fairly  uniform  In  regions  extending  across  a  large 
part  of  the  local  Jet  diameter.  Lastly,  the  amount 
of  unmlxed  ambient  fluid  on  the  jet  centerline  Is 
found  to  vary  roughly  periodically  with  a  period 
approximately  equal  to  the  local  characteristic 
large  scale  time  of  the  riow.  These  results  suggest 
that  large  scale  transport  mechanisms,  displaying  a 
characteristic  organization,  play  an  important  role 
In  entrainment  and  mixing  In  the  far  field  of 
turbulent  Jets. 


Introduction 

Experimental  Investigations  over  the  past 
decade  have  demonstrated  that  entrainment  and  mixing 
In  fully  turbulent  shear  layers  Is  character lzed  by 
organization  that  results  from  the  dynamics  of  large 
scale  vortical  motions.  These  large  scale  motions 
have  been  shown  to  transport  unmlxed  fluid,  from 
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both  free  streams,  across  the  entire  extent  of  the 
layer  (  Konrad'  1976,  Breldenthal*  1978  ). 
Furthermore,  as  a  consequence  of  these  motions,  the 
probability  of  finding  mixed  fluid  within  a  given 
range  of  compositions  Is  essentially  uniform  across 
the  entire  layer  (  Konrad1  1976,  Fiedler*  1974, 
Koochesfahanl '  1989,  Koochesf ahani  &  Dimotakis* 
1984).  On  the  basis  of  these  results  and  many 
others,  It  is  now  generally  recognized  that  models 
of  entrainment  and  mixing  in  turbulent  shear  layers 
probably  need  to  Incorporate  features  of  these  large 
scale  motions. 

Such  large  scale  organized  mixing  may  also  be 
present  In  other  free  turbulent  shear  flows, 
Including  the  axlsymmetrlc  free  Jet.  The  existence 
of  organized  vortical  structure  in  the  near  field  of 
turbulent  Jets  has  been  recognized  for  some  time 
(e.g.  Bradshaw  et  al‘1964,  Crow  &  Champagne7  1971). 
The  prevailing  view,  however,  Is  that  this  near 
field  structure  does  not  survive  In  any  organized 
form  beyond  the  first  few  jet  diameters.  Mixing  in 
the  far  field  of  the  Jet  is  classically  viewed  as  a 
stochastic  process,  Involving  transport  by  small 
scale  eddies  driven  by  gradients  In  the  mean 
profiles,  characteristically  lacking  any  persistent 
large  scale  organization.  This  picture  has  arisen 
largely  from  point  measurements  and  flow 
visualization  techniques  which  are  unable  to  reveal 
the  internal  structure  of  the  mixing. 

However,  the  results  of  several  recent 
Investigations  (  Tso  et  al*  1981,  Dimotakis  et  al* 
1983a,  Dimotakis  et  al,#  1983b,  Dahm  et  al"  1984  ) 
have  suggested  that  the  relevant  transport 
mechanisms  In  the  far  field  of  turbulent  Jets  may  be 
quite  different  from  this  classical  picture  and  may 
exhibit  a  large  scale  organization  similar  in  many 
respects  to  that  found  In  shear  layers.  To 
Investigate  this  prospect  further,  the  measurements 
presented  here  were  conducted  to  study  several 
aspects  of  entrainment  and  mixing  In  turbulent  Jets. 
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Experimental  Methods 

Three  sets  of  measurements  of  entrainment  and 
mixing  in  the  self-similar  far  field  of  an 
axisymmetric  momentum-driven  free  turbulent  jet  in 
water  (  Sc  » 600  )  are  described  below.  These 
measurements  were  all  carried  out  in  a  water 
facility,  documented  previously’ • 1 1 ,  using  laser 
induced  fluorescence  (LIF)  techniques. 


I.  Flame  length  fluctuations  of  reacting  jets 

The  chemically-sensitive  LIF  technique- . 1 0 . 1 1 
was  used  in  conjunction  with  an  (isothermal) 
acid-base  reaction  to  determine  length  and  time 
scales  for  the  flame  length  fluctuations  of 
chemically  reacting  jets. 

Briefly,  an  alkaline  jet  premixed  with  a  small 
amount  of  a  pH-sensitive  laser-fluorescent  dye  was 
discharged  through  an  axisymmetric  nozzle  into  a 
large  acidic  reservoir.  By  appropriately  selecting 
the  initial  absolute  concentrations  of  the  acid  and 
base  solutions,  the  pH-sensltive  characteristics  of 
the  dye  could  be  exploited  so  that  mixed  fluid 
illuminated  by  the  laser  excitation  would  fluoresce 
only  if  the  local  mass  ratio  of  ambient  fluid  to  jet 
fluid  mixed  at_  the  molecular  scale  did  not  exceed  a 
threshold  value,  referred  to  here  as  the  reaction 
stoichiometric  ratio  $.  This  ratio  could  be 
selected  by  choosing  the  relative  initial 
concentrations  of  the  acid  and  base  solutions.  The 
fluorescence  transition  was  reversible  and  occured 
over  a  very  narrow  range  of  the  mixture  ratio.  By 
arranging  the  laser  light  to  form  either  a  single 
thin  sheet  containing  the  jet  axis  or  a  pair  of 
mutually  orthogonal  sheets  both  containing  the  Jet 
axis,  the  jet  fluid  in  the  plane  of  the  laser 


sheet(s)  that  had  not  yet  mixed,  on  a  molecular 
scale,  with  ambient  fluid  in  excess  of  the  selected 
stoichiometric  ratio  could  be  visualized  directly. 

Earlier  measurements  (Dahm  et  al11  198*4)  using 
this  technique  have  demonstrated  a  nearly  periodic- 
fluctuation  of  the  "flame"  length  of  such  reacting 
jets.  A  composite  sequence  showing  these 
fluctuations,  consisting  of  every  fifth  frame  from 
an  excerpt  of  a  high  speed  motion  film  of  such  a 
chemically  reacting  turbulent  jet  visualized  using 
this  technique,  is  shown  in  figure  1.  In  each 
frame,  the  fluorescence  from  each  of  a  pair  of 
mutually  orthogonal  laser  sheets,  both  containing 
the  jet  axis,  was  imaged  by  a  system  of  mirrors 
without  perspective  distortion  onto  the  film  plane 
to  allow  a  simultaneous  view  of  mixing  in  the  two 
thin,  orthogonal  slices  through  the  jet. 

While  the  effects  of  heat  release  and  density 
changes  that  accompany  combustion  are  not  addressed 
in  experiments  with  such  simple  acid-base  reactions 
in  water,  experience  in  turbulent  mixing  layers  has 
shown  that  the  basic  ideas  suggested  by  such 
investigations  remain  applicable  when  combustion  is 
present.  Indeed,  when  properly  normalized,  the  mean 
flame  lengths  of  momentum-driven  turbulent  jets  in 
water11!12  and  in  highly  exothermic  gas  flames12-17 
as  well  as  in  liquid-vapor  condensation  systems18-20 
are  in  good  agreement,  as  shown  in  figure  2.  Such  a 
correlation  was  demonstrated  by  Avery  &  Faeth18 
(1971*).  Broadwell21  (1982)  uses  a  similar  form,  for 
water  and  gas  flames,  to  argue  that  the  rate  of 
molecular  mixing  is  set  by  the  rate  of  entrainment 
for  the  jet.  Figure  2  then  demonstrates  that,  at 
Reynolds  numbers  sufficiently  large  for  mixing  to 
become  entrainment-limited,  the  flame  length  of  any 
reacting  jet  is  determined  solely  from  the  generic 
scaling  laws  for  jet  entrainment  expressed  in 
momentum  variables. 


Figure  1.  Composite  time  sequence  from  a  motion  film  of  a  chemically  reacting  turbulent  jet, 
visualized  using  the  dual -orthogonal  planar  LIF  technique,  showing  the  flame  length 
fluctuations;  ^  *  18,  Re  -  10,000. 


The  correlation  In  figure  2  arises  from  the 
normalization  with  d»,  the  jet  "momentum  diameter", 
defined  as 


where  m0  is  the  jet  exit  mass  flux,  J0  Is  the  Jet 
momentum  flux,  and  9.  is  the  ambient  fluid  density. 
This  momentum  diameter  was  Introduced  In  a  limited 
context  by  Thrlng  4  Newby1 2  (1952),  and  has  been 
used  by  Becker  et  al“.“  (1978)  In  gas  flames. 
However,  figure  2  suggests  a  general  scaling  with  d* 
for  all  momentum  Jets,  and  Indicates  a  proper 
dimensionless  form  for  the  axial  coordinate  x  of 
the  far  field  of  all  momentum-driven  turbulent  jets 


In  the  present  experiment,  the  Instantaneous 
turbulent  flame  length  was  measured  directly,  using 
a  coordinate  digitizer,  from  each  of  over  2900 
Individually  projected  film  frames  of  a  high  speed 
motion  film  of  such  a  reacting  jet  at  Re  -  10,000  and 
♦  -15.  The  resulting  data  spanned  roughly  50 
periods  of  the  flame  fluctuation.  Length  and  time 
scales  for  these  fluctuations  were  determined  from 
the  measurements  and  compared  with  the  local 
characteristic  large  scale  length  and  time  of  the 
flow,  respectively  the  local  Jet  visual  diameter 
D  -  0 . 9  x  at  the  mean  flame  tip  location  x,  and  D/Uc 
where  U  .  6  (J„/m  )  v  1  Is  the  mean  centerline 
velocity  at  the  flame  tip. 
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Figure  2.  Dimensionless  mean  length,  L,  of 
momentum-driven  turbulent  Jet  names  as  a  function 
of  their  reaction  stoichiometric  ratio,  $,  for 
Isothermal  reactions  In  water,  exothermic  gas 
reactions  and  liquid-vapor  condensation  reactions. 
(Buoyancy  effects  are  present  to  various  degrees  In 
some  of  these  data.) 


II.  Radial  profiles  of  concentration 
in  non-reacting  Jets 

Instantaneous  radial  profiles  of  Jet  fluid 
concentration  In  a  non-reacting  Jet  were  measured 
with  high  resolution  using  a  passive  LIF  Imaging 
technique  In  conjunction  with  high  speed  digital 
data  acquisition  and  digital  Image  display  methods. 
The  technique  was  similar  to  earlier  measurements  of 
concentration  In  the  plane  shear  layer 
(Koochesfahanl  4  Dimotakis’  1989). 


the  fluorescence  Intensity.  A  green  filter 
effectively  eliminated  directly  scattered  laser 
light  from  any  particulates  In  the  flow. 

Radial  profiles  of  fluorescence  Intensity  were 
measured  at  p  300  for  Re  -  1500  and  Re  -  5000,  on 
either  side  of  the  Re »  3000  asymptotic  threshold 
(  Dahm  et  al 1 1  1989  )  at  which  mixing  in  the  far 
field  of  such  Jets  In  water  becomes 
entrainment-llmlted  and,  therefore,  the  flame  length 
becomes  Independent  of  Reynolds  number. 


The  laser  fluorescent  dye  was  again  premixed 
with  the  Jet  fluid  (water),  and  discharged  Into  the 
ambient  reservoir  fluid  (water).  In  this  case,  no 
chemical  reaction  was  Involved.  The  beam  (519.5  nm) 
from  an  argon  ion  laser  was  collimated  and  oriented 
to  cross  the  Jet  radially,  perpendicular  to  the  Jet 
axis.  The  fluorescence  from  dye  along  the  beam  was 
Imaged  onto  a  1029-element  Retlcon  self-scanning 
linear  photodiode  array,  the  sequential  pixel  output 
of  which  gave  consecutive  instantaneous  profiles  of 


The  photodiode  array  was  clocked  at  a  pixel 
rate  of  256  kHz,  corresponding  to  a  scan  rate  of  235 
scans/sec  and  a  scan  time  of  9.3  msec,  which  was 
several  orders  of  magnitude  faster  than  any  relevant 
local  time  scales  in  the  flow  at  either  Reynolds 
number.  Consequently,  temporal  resolution  was  not 
an  Issue  In  these  measurements. 

Each  of  the  1029  photodiode  array  elements 
(pixels)  had  a  sensitive  detector  aperture  of 
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25uin  «  26um.  The  optics  used  produced  an  image  ratio 
of  26:1,  and  the  beam  diameter  varied  from  roughly 
0.5mm  to  about  1mm,  thereby  defining  the  measurement 
volume  for  each  pixel.  The  resulting  resolution  was 
sufficient  to  resolve  even  the  smallest  turbulent 
scales  of  motion  at  the  lower  Reynolds  number.  At 
neither  Reynolds  number  could  the  strain-limited 
diffusion  layer  thickness  (  Batchelor  scale  )  be 


-1/2 

resolved,  which  is  of  the  order  Sc  smaller  than 


the  Kolmogorov  scale. 


The  array  output  was  digitized  through  a  single 
channel  high-speed  8-bit  A/D  converter  and  recorded 
on  a  computer  disk.  Since  the  scan  time  was  fast 
relative  to  the  local  time  scales  in  the  flow,  and 
in  order  to  keep  the  amount  of  data  collected 
manageable,  only  every  16th  scan  of  the  array  was 
actually  recorded  on  the  disk  at  Re*  1500,  and  every 
6th  scan  at  Re*  5000.  Each  measurement  consisted  of 
4000  recorded  scans  of  the  array  (in  excess  of  4 
million  measurements)  spanning  approximately  12 
local  characteristic  large  scale  flow  times,  this 
limit  being  Imposed  by  the  local  flow  duration. 


The  measured  pattern  of  dark  noise  and 
background  illumination  was  subtracted  from  these 
data.  The  resulting  measured  fluorescence  intensity 
profiles  could  then  be  converted  to  profiles  of  Jet 
fluid  concentration  a3  described  below  (  see  also 
Koochesfahani  4  Dimotakis5  1964). 


The  intensity  of  fluorescence  measured  by  the 
photodiode  array,  F(E,t),  was  locally  related  to  the 
Instantaneous  dye  concentration  and  beam  Intensity 
at  every  point  along  the  beam  as 


F(E.t)  -  h(£)  •  J [c( E. t )  J  1(5,  t) 


where  Elsa  coordinate  along  the  beam  in  the 
direction  of  propagation  and  c(E.t)  is  the 
Instantaneous  dye  concentration  profile  along  the 
beam.  Here,  I(E,t)  is  the  instantaneous  beam 
intensity  profile,  j[c]  relates  the  local  intensity 
of  fluorescence  to  the  local  beam  intensity,  and 
h(E)  is  the  transfer  function  relating  the 
fluorescence  intensity  to  that  measured  by  the 
photodiode  array. 


The  instantaneous  beam  intensity  profile  I(E,t) 
was,  in  turn,  related  to  the  dye  concentration 
profile  c(E.t)  through  the  attenuation  integral 


I(E.t) 


>(  -  j[cU,t)  ]  dA  ) 


where  IQ(t)  is  the  beam  intensity  at  E-0  and  A(E) 
is  the  beam  cross  sectional  area  profile.  Over  the 
range  of  concentrations  in  this  experiment,  j[c]  was 
determined  by  a  calibration  to  be  linearly  related 
to  the  dye  concentration  as 


>  v;.-  v 


j[c]  *  C  •  c 


where  the  constant  c.  related  to  the  scattering 
cross  section  of  the  dye,  was  determined  from  the 


calibration. 


Defining  an  overall  transfer  function  H(E)  to 
include  the  beam  area  variation,  namely 


H(E)  -  h( E) 


gave  the  instantaneous  fluorescence  intensity 
profile  measured  by  the  photodiode  array  as 


F(S,t)  -  IQ(t)  H(E)  C c(E.t)  exp(  c(A,t)dAj. 


The  transfer  function  H( E)  was  determined  from  an 
in  situ  measurement  of  F ( E )  at  a  uniform  dye 


concentration  c1  as 


H ( E )  /  H( E* ) . 


F(E>  '  Ft E* ) ]  exp  (  t  c  (e*E*)  ) 


where  E*  “as  an  arbitrary  reference  location.  H(E) 
represents,  therefore,  an  accounting  of  the 
collective  non-idealities  in  the  optical  system  and 
the  array. 


Normalizing  each  instantaneous  measured 
fluorescence  intensity  profile  by  it3  value  at  E*. 
where  the  dye  concentration  was  at  a  constant 
reference  value  c*,  gave  the  dye  concentration 
profile  relative  to  this  reference  as 


c(E.t) 


F(E.t) /  F(E*,t) . 


[  H( E)  /  H ( E* )  J 
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A  small  vessel,  positioned  about  one-and-a-half 
local  jet  diameters  from  the  centerline,  containing 
a  circulating  dye  solution  at  a  known  concentration, 
was  also  imaged  onto  the  array  to  set  this  reference 
condition. 


Denoting  the  known  concentration  of  dye  issuing 


from  the  Jet  exit  as  cq  defined  the  Instantaneous 


jet  fluid  concentration  profile  c.{£,t)  in  terms  of 
the  dye  concentration  profile  as 


CjCE.t) 


c(E.t) 


The  Jet  fluid  concentration  profile  along  the  beam 
corresponding  to  each  measured  fluorescence 
intensity  profile  could  then  be  determined  as 
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The  resulting  instantaneous  profiles  of  jet 
fluid  concentration  were  expressed  in  the 
conventional  radial  self-similar  form  g(n,t),  where 
n*r/x,  by  normalizing  each  individual  profile  with 
the  mean  centerline  concentration  as 

c.(n,t) 

g(n,t)  -  j-  —  . 

Cj(0) 


and  $-15.  The  nearly  regular  length  fluctuations 
of  figure  1  are  evident  in  these  data.  Defining 
(albeit,  somewhat  subjectively)  fluctuation  events, 
delimited  by  tic  marks  in  figure  3,  allows  a  length 
and  time  scale  to  be  determined  for  each  fluctuation 
event . 

Histograms  of  the  resulting  fluctuation  length 
and  time  scales,  normalized  by  the  corresponding 
local  characteristic  large  scale  quantities,  a  re 
shown  in  figure  The  histograms  indicate  tnat  the 
flame  length  fluctuations  of  reacting  jeta  occur  on 
a  length  and  time  scale  approximately  equal  to  the 
local  characteristic  large  scales  of  the  flow. 


III.  Unmixed  ambient  fluid  on  the  jet  axis 

The  probability  of  unmixed  ambient  fluid 
appearing  at  a  point  on  the  Jet  axis  was  measured  as 
a  function  of  time  for  Re-  1500  and  Re  -  5000  using  a 
passive  LIF  technique. 

The  laser  beam  was  oriented  coincident  with  the 
jet  centerline.  The  laser  dye  was  again  premixed 
with  the  Jet  fluid,  and  the  fluorescence  from  a 
segment  along  the  beam  centered  at  y  -  300  and 
extending  from  285  <  x  <  3'5  (roughly  one-fifth  the 
3cale  of  the  local  Jet  diameter)  was  imaged  onto  the 
1024-element  Reticon  linear  photodiode  array.  A 
green  filter  was  again  used  to  block  directly 
scattered  laser  light.  As  before,  the  array  was 
scanned  at  256  kHz  and  the  output  from  every  16th 
scan  at  Re  -  1500  and  every  6th  scan  at  Re  -  5000  was 
recorded.  The  data  for  each  of  these  measurements 
consisted  or  tt000  scans  of  the  array,  again  spanning 
approximately  12  local  characteristic  large  scale 
flow  times. 

In  this  case,  however,  it  was  not  necessary  to 
convert  the  measured  fluorescence  intensities  to  jet 
fluid  concentration,  as  was  done  in  the  previous 
experiment,  but  only  to  distinguish  between  mixed 
and  unmixed  fluid.  Consequently,  a  reference 
condition  was  not  used  in  this  measurement.  The 
threshold  between  mixed  and  unmlxed  fluid  was 
determined  from  the  peak  dark  noise  in  the  array. 
The  fraction  of  the  1024  pixels  Imaging  unmixed 
ambient  fluid,  corresponding  to  the  probability  of 
finding  unmlxed  fluid,  could  subsequently  be 
determined  as  a  function  of  time. 


Results  and  Discussion 


These  fluctuations  can  be  interpreted  in  the 
context  of  the  instantaneous  concentration  field  of 
the  jet.  Their  length  and  time  scaling  show  that 
the  mixed  fluid  composition  in  a  large  region  of  the 
Jet  crosses  the  reaction  stoichiometric  ratio  within 
one  local  cnaracteristic  large  scale  time, 
suggesting  that  this  fluid  is  within  a  fairly  narrow 
range  of  compositions.  Such  regions  of 
qualitatively  uniform  jet  fluid  concentration  have 
been  inferred  from  still  photographs  of  planar 
fluorescence  in  jets  (Dimotakis  et  al 1 0  1933b).  The 
flame  length  fluctuations  are,  therefore,  directly 
indicative  of  a  nearly  periodic  large  scale 
organization  of  mixing  in  the  jet  far  field. 
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I.  Flame  length  fluctuations  of  reacting  jets 

Figure  3.  Instantaneous  turbulent  flame  length 
vs  .  time;  >  -  1 1 ,  Re  -  1 0,000 . 


The  instantaneous  turbulent  flame  length  as  a 
function  of  t ,  me  Is  shown  in  figure  3  for  Re  -  10,000 
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concentration  (in  self-similar  form)  compared  with 
other  experiments. 


Figure  M,  Histograms  of  length  and  time  scales  for 
the  flame  length  fluctuations  in  figure  3. 
normalized  by  the  corresponding  local  character ist i c 
large  scale  flow  quantities. 


Hussain23  (1983)  has  suggested  that  the 
significance  of  organization  in  turhulent  shear 
flows  ha3  been  overemphasized.  The  flame  length 
fluctuations,  however,  suggest  that  this  large  scale 
organization  plays  a  crucial  role  in  the 
entrainment,  mixing  and  chemical  reaction  processes 
in  the  jet.  It  is  noteworthy  that  similar  flame 
length  fluctuations  have  also  been  reported  in 
combustion  studies  in  buoyant  plume  flows 
(ZuKOSki  et  al3'  198M). 


II.  Radial  profiles  of  concentration 
in  non-reacting  Jets 

The  mean  radial  concentration  profile  at 
Re -  5000,  determined  by  an  ensemble  average  of  the 
individual  instantaneous  profiles  measured  in  these 
experiments.  Is  compared  In  figure  5a  In  the  radial 
similarity  variable  n-r/x  with  the  profiles 
reported  In  other  experiments33'31.  Note  that  this 
form  of  the  radial  variable  n  does  not  force  the 
profiles  to  match  at  the  Jet  half-width,  and  rhus 
provides  a  more  rigorous  comparison  than  does 
normalization  of  the  radial  coordinate  by  the 


Figure  5b.  Radial  profile  of  rms  concentration 
fluctuations  (in  self-similar  form)  compared  with 
other  experiments. 


concentration  half-width.  Similarly,  the  radial 
profile  of  rms  concentration  fluctuations  is  shown 
in  figure  5b.  However,  the  principal  emphasis  in 
the  present  experiments  is  on  high  resolution 
measurements  at  a  large  number  of  points  across  the 
Jet  and,  consequently,  the  duration  over  which  these 
measurements  extend  is  relatively  short. 
Nevertheless,  despite  the  limited  statistical 
significance  of  mean  profiles  obtained  from  these 
data,  the  comparisons  in  figures  5a  and  5b  are  quite 
good  and  serve  to  validate  the  measurements. 

Several  representative  instantaneous  radial 
concentration  profiles  at  Re -  5000,  each  covering 
the  entire  radial  extent  of  the  jet,  are  shown  in 
figures  6a-d. 

Perhaps  the  most  significant  feature  of  the 
Instantaneous  profiles  13  how  little  they  resemble 
the  mean  profile.  Most  of  the  instantaneous 
profiles  are  either  of  a  "top-hat"  type,  such  as 
figures  6c  and  6d,  where  the  mixed  fluid  across  the 
entire  jet  diameter  is  within  a  fairly  narrow  range 
of  concentrations  (relative  to  what  the  mean  profile 
suggests),  or  of  a  "two-level"  type,  such  as 
figures  6a  and  6t,  where  the  mixed  fluid  in  a  large 
region  near  the  jet  axis  is  within  one  such  narrow 
range  of  concentrations,  while  the  mixed  fluid 
across  the  remainder  of  tne  jet  is  within  a 
second,  lower  range  of  concentrations. 
Broadwell  4  Bre identhal 1 2  (1982)  nave  noted  that 
such  qualitative  uniformity  in  the  mixed  fluid 
composition  can,  in  retrospect,  also  be  recognized 
in  concentration  profiles  from  plane  turbulent  jets 
by  Uberol  4  Singh11  (1975),  but  that  these  authors 
do  not  remark  about  the  shape  of  their  profiles. 

Furthermore,  the  profiles  in  figure  6  show  that 
unmixed  ambient  fluid  can  be  found  throughout  the 
jet.  This  has  been  observed  in  still  photographs  of 
planar  fluorescence  in  jets  (Dimotakis  et  al* 
1983a),  and  supports  the  proposal  that  the  classical 
notion  of  entrainment  and  mixing  by  stochastic  small 
3cale  eddies  may  not  be  a  good  picture  of  the 
relevant  transport  mechanisms  in  the  Jet  far  field. 

The  mean  concentration  profile  of  figure  5a 
then  results  from  an  average  of  the  mixed  fluid 
profiles  described  above  and  the  unmixed  fluid 
within  the  jet,  but  is  itself  a  poor  representat ion 
of  the  compositions  within  the  Jet.  Thi3  contrasts 
with  the  classical  picture  of  turbulent  mixing,  in 
which  instantaneous  concentration  profiles  are 
presumed  to  qualitatively  resemble  the  mean  profile 
with  small  turbulent  fluctuations  about  the  mean, 
and  may  have  particularly  important  implications  for 
the  modeling  and  prediction  of  chemical  reactions 
between  the  mixed  and  unmixed  fluids. 
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Figure  6.  Representative  instantaneous  radial 
profiles  of  concentration  across  the  jet;  x"  300, 
Re  -  5000. 


1  ^  k'lVi  LT  A-  A  V  1 


The  instantaneous  concentration  profile  shapes 
in  figures  6a-d  as  well  as  the  observed  flame  length 
fluctuations  of  figures  2  and  3  can  be  interpreted 
in  terms  of  the  simple  conceptual  picture  of  mixing 
in  the  jet  far  field  shown  in  figure  7,  described  by 
Dimetukis  et  al10  (1933b).  In  this  picture,  large 
regions  within  the  jet  contain  mixed  fluid  in  a 
fairly  narrow  range  of  compositions,  intertwined 
with  unmixed  ambient  fluid.  As  a  result, 
concentration  profiles  across  the  jet  will  nave 
either  the  "top-nat"  or  ".wo-level"  shape  noted 
above,  while  the  mixed  fluid  concentration  along  the 
axial  direction  decreases  in  a  roughly  step-wise 
fashion,  following  the  reruir.it.-  x_1  form  in  the 
mean,  and  giving  rise  to  the  flame  length 
fluctuations  in  a  reacting  jet. 


A  method  for  simul  t  -one  uusly  interpreting  many 
of  these  measured  instant  arn-ous  concentration 
profiles  is  shown  in  figure  5.  These  pictures  are 
similar  to  earlier  displays  of  LIF  data  in  a  plane 
mixing  layer  (Koochesf  uhani  \  Oiir.otikin4  193H. )  They 
are  photographs  of  a  digital  image  display  screen, 
on  which  trie  digitized  concentration  data  from  the 
512  array  elements  centered  on  the  jet  axis  were 
displayed  for  500  consecutive  recorded  array  scans. 


Figure  7.  Idealized  conceptual  picture  of  the 
instantaneous  concentration  field  of  the  turbulent 
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Figure  3.  r-t.  reeon:;?  ruet  ion  o!  ine  turbulent  j 
coneerKr-i i  i on  fieiO  frora  000 
eon-.-'fit r.if.  ion  [’rofilrn;  x  =  W0. 


The  256  different  grey  levels  are  assigned  to  denote 
the  various  coneentru t ion  levels,  with  pare  while 
representing  g(n,t)  =1.5  and  decreasing  linearly  to 
pure  bl3ck  denoting  unmixed  ambient  fluid,  namely 
g(n,t)=0.  These  are,  therefore,  r-t  pictures  of 
the  concentration  along  a  line  across  the  jet.  It 
must  be  emphasised  tint  these  are  not  equivalent  to 
still  phologr  iph.s  of  the  concentration  field  in  a 
plane  intersecting  the  jet,  as  becomes  evident  near 
the  jet  edges. 

It  is  apparent  in  the  data  in  figure  S  that,  at 
both  Reynolds  numbers,  unmixed  ambient  fluid  is 
transported  throughout  the  entire  radial  extent  of 
the  jet.  This  figure  also  suggests  that  the 
topology  of  a  stoichiometric  surface  between  the 
mixed  and  unmixed  fluids,  commonly  used  in 
combustion  studies,  may  be  considerably  more  complex 
than  that  which  might,  be  suggested  by  a  conventional 
small  scale  trar-pcrt  picture  of  jet  mixing. 
Furthermore,  as  the  individual  concentration 
profiles  in  figure  6  suggested,  the  data  in  figure  8 
show  that  the  instantaneous  mixed  fluid  composition 
can  be  fairly  uniform,  relative  to  what  the  mean 
concentration  profile  would  suggest,  across  a  large 
part  of  the  local  jet  diameter. 

A  composite  sequence  from  consecutive 
photographs  of  the  typo  displayed  in  figure  S  is 
shown  for  each  Reynolds  number  in  figure  9.  Each  of 
these  sequences  displays  roughly  2  million 
individual  measurements.  This  method  of  displaying 
tne  nciairuJ  profiles  allows  a  direct  interpretation 
of  tne  tine-varying  concentration  field  in  these 
data.  Tine  transport  of  unmixed  ambient  fluid 
throughout  the  jet,  as  well  as  regions  of  fairly 
uniform  mixed  fluid  concentration  extending  across  a 
large  part  of  tile  local  jet  diameter,  can  again  be 
discerned  in  these  pictures.  Note  that  a  length 
increment.  along  the  jet  centerline  in  the 
pseudo-axial  (temporal)  direction  can  be  made 
roughly  equivalent  to  a  corresponding  increment  in 
tin1  radial  direction  by  a  50%  stretching  of  tne 
temper  :  1  d i ree t i on . 


III.  Unmixed  ambient,  fluid  on  the  jet  axis 


The  fraction  of  tin-  array  detecting  unmixed 
ambient  fluid  along  the  imu.j-.J  segment  of  the  jet 
axis  Is  shewn  as  a  function  of  tint-*  for  both 


Reynolds  numbers  in  figure  to. 


These  data  show  that  the  probability  of  finding 
utv'.i  ambient  flui  d  at  a  point  on  the  jet  axis 
show:,  a  m  irked  increase  at  roughly  period!  : 
intervals.  Furthermore,  as  indicated  by  the  time 
scale  accompany i ng  each  trace  these  increases  are 
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Fig.  10a  Re  -  1500. 


scale,  undergoes  a  shift  toward  higher  wavenumbers 
as  the  Reynolds  number  is  increased,  and  that 
consequently  less  unmlxed  fluid  is  detected.  In 
other  words,  the  flow  at  Re*  1500  may,  in  this 
respect  at  least,  be  qualitatively  different  than 
that  at  Re  *  5000,  an  idea  that  is  supported  by  the 
observation  that  reacting  jets  in  water  attain  their 
asymptotic  flame  lengths  at  Re»  300011.  This 
question,  however,  must  ultimately  be  resolved  using 
other  techniques. 

Conclusions 


Po  1 


The  results  of  these  measurements,  together 
with  previous  investigations’ .  1  ° , 1  1 ,  indicate  that 
large  scale  transport  mechanisms  play  an  important 
role  in  entrainment  and  mixing  in  the  far  field  of 
turbulent  jets,  and  that  the  jet  far  field  is 
characterized  by  large  scale  organization. 
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Fig.  10b  Re  =  5000. 


Figure  10.  Fraction  of  the  imaged  segment  on  the 
jet  centerline  detecting  unmixed  ambient 
fluid  vs.  time;  x  -  300. 


typically  separated  by  about  one  local 
characteristic  large  scale  flow  time.  This  roughly 
periodic  Increase  in  the  probability  of  finding 
unmixed  ambient  fluid  on  the  jet  axis  is  consistent 
with  the  large  scale  entrainment  mechanism  proposed 
by  Dimotakis  et  al’  (1933a)  and  with  the  conceptual 
picture  of  jet  mixing  in  figure  7. 

In  figure  10,  the  time-averaged  probability  of 
detecting  unmlxed  ambient  fluid  at  a  point  on  the 
jet  axl3  is  29%  at  Re-  1500,  while  at  Re  -  5000  this 
value  is  3%.  This  difference  may,  to  a  certain 
extent,  be  ascribed  to  the  decrease  in  relative 
resolution  of  the  measurement  at  the  higher  Reynolds 
number.  The  effect  of  resolution  on  these 
measurements,  however,  may  be  determined  by 
averaging  the  LIF  intensity  over  two  or  more  pixels 
before  applying  the  mlxed/unmixed  criterion.  The 
result  of  such  a  procedure  indicates  that,  when  the 
effective  resolution  at  Re -1500  is  made  equivalent 
to  that  in  the  original  measurement  at  Re-  5000,  the 
time  averaged  probability  of  detecting  unmixed  fluid 
is  still  24%.  While  we  might  be  led  to  conclude 
from  this  that  less  unmlxed  ambient  fluid  reaches 
the  centerline  at  the  higher  Reynolds  number,  it  is 
also  possible  that  the  spectrum  of  scales  at  which 
unmixed  fluid  is  found,  relative  to  the  Kolomogorov 


Instantaneous  profiles  of  concentration  across 
the  jet  show  that,  at  any  instant,  the  mixed  fluid 
concentration  can  be  fairly  uniform  within  large 
regions  in  the  jet,  and  that  unmixed  ambient  fluid 
can  be  found  throughout  the  jet.  These 
characteristics  of  the  instantaneous  profiles  as 
well  as  the  measurements  of  unmixed  fluid  within  the 
jet  are  consistent  with  a  3imple  conceptual  picture 
of  large  scale  organization  of  entrainment  and 
mixing  in  jets.  For  the  case  of  chemically  reacting 
jets,  this  organization  manifests  itself  as  a  nearly 
periodic  large  scale  fluctuation  of  tne  flame 
length,  indicating  the  important  role  that  such 
large  scale  organization  plays  in  nixing  in  the  far 
field  of  turbulent  jets. 
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